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Basic program linearizes thermistors 



Mick Murray 

Lab-Line Instruments Inc, Melrose Park, IL 

You can reduce the nonlinearity of a negative-tempera- 
ture-coefficient (NTC) thermistor (TH-1) by adding a 
series resistor (RS) and a parallel resistor (R). Fig 1 
contains the resulting network, designated R E q. To 
define optimum values for the resistors, the Basic 
program of Listing 1 uses an iterative process, select- 
ing from standard values in data statements 480 
through 560. 

Using a 3-point straight-line approach, the program 
can linearize any NTC thermistor. First, the program 
prompts you for the lowest and highest temperatures of 
interest and the corresponding thermistor resistance 



(which you can find in the manufacturer's data sheet). 
Next, it calculates the exact mid-range temperature 
and prompts you for the resistance at that point as well. 
You divide each resistance value by lk before entering 
it in the program. 

In the example of Fig 1, thermistor TH-1 varies 127 
to 2.64 kft over the temperature range of 20 to 130°C. 
Resistance change in the lower half of this range (20 to 
75°C) is 9.3 times greater than that in the upper half. 
By comparison, the computer-selected network R K q 
varies only 1062 to 977fl over the same temperature 
range, but the change in each half of the range is nearly 
equal. (Note, however, that the network's linearity 
error increases as you widen the specified temperature 
range.) 

Text continued on pg 222 



Reo 



1.03k 

S 

g 1.02k 




TEMPERATURE ("CI 



Fig 1 — The Ilasic program in Listing 1 selects the standard values for resistors R and RS that provide the most nearly linear relationship 
between Rro and temperature. 
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LISTING 1 — LINEARIZATION PROGRAM 



10 REM ******************************** 

20 REM 

30 REM THERMISTOR LINEARIZATION PROGRAM 
40 REM COPYRIGHT 1986 MICK MURRAY 
50 REM LAB-LINE INSTRUMENTS, INC. 
60 REM MELROSE PARK, IL. 60160 
70 REM 

80 REM ******************************** 

90 A$=STRING$( 70, "*" ) 

100 FOR T=l TO 24:PRINT:NEXT 

110 INPUT"INPUT LOW TEMP ERATURE : " ; L 

120 INPUT"INPUT LOW TEMP RE S I ST ANCE : " ; RL 

130 INPUT"INPUT HIGH TEMPERATURE :"; H 

140 INPUT"INPUT HIGH TEMP RE S I ST ANC E : " ; RH 

150 M=( (H-L)/2)+L 

160 PRINT "THE CALCULATED MID-POINT TEMPERATURE = " ; M ; "DEGREES. " 

170 PRINT"INPUT THERMISTOR RESISTANCE AT ";M;" DEGREES:"; 

180 INPUT RM:ST=1000 

190 REM DE=DELTA, ST=STORE 

200 PRINT A$:PRINT:PRINT 

210 DIM A(96) :DIM B(96) 

220 FOR X=l TO 96 : READ A( X ) : B ( X ) = A( X ) : NEXT X 
230 FOR 1=1 TO 96 
240 FOR J = l TO 96 

2 50 L0= ( ( RL+A( I ) )*B( J) ) / ( RL + A( I ) +B ( ) ) 

2 60 MI = ( (RM+A( I ) )*B( J) )/(RM+A( I ) + B( ) ) 

2 70 HI = ( (RH+A( I) )*Bf J) ) / ( RH+A( I ) +B ( J) ) 

280 X=ABS(HI-MI) 
290 X1=ABS(MI-L0) 
300 DE=ABS(X1-X) 

310 IF DE<=ST THEN ST - DE : RS = A ( I ) : R= B ( J ) : GOSIIB 380 

320 NEXT J 
330 NEXT I 

340 IF SD = THEN GOTO 580 

350 PRINT" END OF CALCULATIONS." 

360 PRINT CHR$( 7 ) ; : F0RP=1 T01 000 : NEXT:G0T036O 

370 END 

380 PRINT "BEST SO FAR:" 

390 PRINT"RS = " ; A( I ) ; " OHMS , AND R = " ; B ( ) ; "OHMS . " 
400 PRINT" *** ( AT" ; L ; " DEGRE ES , Rth = " ; LO ; "OHMS . ) " 
410 PRINT" *** (AT" ;M; "DEGREES, Rt.h = " ; M I ; "OHMS . ) " 
420 PRINT" *** (AT" ;H; "DEGREES, Rth = " : H I ; "OHMS . ) " 
430 PR I NT" DELTA1 =";X;"OHMS, AND DELTA2 = " ; X 1 ; "OHMS . " 
440 HL=ABS(HI-LO) 

450 PRINT "THE TOTAL CHANGE IN RESISTANCE FROM " ;L;" TO ";H;" DEGREES I S " ; HL ; "OHMS . " 
460 PR I NT : PR I NT AS : PR I NT"Work i nq . . . " : PR I NT 
470 RETURN 

480 DATA 1.00,1.02,1.05,1.07,1.10,1.13,1.15,1.18,1.21,1.24 

490 DATA 1.27,1.30,1.33,1.37,1.40,1.43,1.47,1.50,1.54,1.58 

500 DATA 1.62,1.65,1.69,1.74,1.78,1.82,1.87,1.91,1.96,2.00 

510 DATA 2.05,2.10,2.15,2.21,2.26,2.32,2.37,2.43,2.49,2.55 

520 DATA 2.61,2.67,2.74,2.80,2.87,2.94,3.01,3.09,3.16,3.24,3.32 

5 30 DATA 3.40,3.48,3.57,3.65,3.74,3.83,3.92,4.02,4.12,4.22,4.32,4.42 

5 40 DATA 4. 5 3, 4. 64,4. 75, 4. 87, 4. 99, 5. 11, 5. 23, 5. 36, 5. 49,5.62,5.76,5, -9 

5 50 DATA 6.04,6.19,6.34,6.49,6.65,6.81,6.98,7.15,7.32,7.5,7.68,7.87 

560 DATA 8.06,8.25,8.45,8.66,8.87,9.09,9.31,9.53,9.76 

570 REM THE FOLLOWING IS FOR THE SECOND DECADE OF RS VALUES: 

580 SD = 1 : RESTORE 

590 FOR N=l TO 96 

600 READ A(N) :B(N)=A(N) :A(N)=10*A(N) 
610 NEXT N:G0T0230 
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With each iteration, the program checks the network 
resistance 'Req at the low, middle, and high tempera- 
tures. It seeks to achieve an equal resistance change in 
the upper and lower halves of the temperature range by 
trying all combinations of R and RS by substituting 
values from the data statements. It prints a result, 
searches for a better result, prints that, and continues 



(Fig 2). Remember to denormalize the final result by 
multiplying the resistance values by lk. Typical execu- 
tion time is about 10 minutes. EDM 



To Vote For This Design, Circle No 748 



INPUT LOW TEMPERATURE: 20 

INPUT LOW TEMP RESISTANCE: 127 

INPUT HIGH TEMPERATURE: 130 

INPUT HIGH TEMPERATURE RESISTANCE: 2.64 

THE CALCULATED MID-POINT TEMPERATURE = 75 DEGREES. 

INPUT THERMISTOR RESISTANCE AT 75 DEGREES: 12.96 

******************************************* 

BEST SO FAR: 

RS = 1 OHMS, AND R = 1 OHMS. 

* * * (AT 20 DEGREES, R EQ = .9922482 OHMS.) 

* * * (AT 75 DEGREES, R EQ = .9331552 OHMS.) 

* * * (AT 130 DEGREES, R EQ = .7844828 OHMS.) 
DELTA1 = .1486723 OHMS, AND DELTA2= .059093 OHMS. 

THE TOTAL CHANGE IN RESISTANCE FROM 20 TO 130 DEGREES IS .2077653 OHMS. 
******************************************* 



(12 ITERATIONS) 



******************************************* 

Working. . . 

BEST SO FAR: 

RS = 8.66 OHMS, AND R = 1.07 OHMS. 

* * * (AT 20 DEGREES, R e q = 1.061627 OHMS.) 

* * * (AT 75 DEGREES, R E0 = 1.019542 OHMS.) 

* * * (AT 130 DEGREES, R EQ = .9774454 OHMS.) 
DELTA1 = .0420962 OHMS, AND DELTA2= .04208481 OHMS. 

THE TOTAL CHANGE IN RESISTANCE FROM 20 TO 130 DEGREES IS .08418102 OHMS. 

******************************************* 

Working . . . 

END OF CALCULATIONS. 



J% 2— This example of Listing I's output shows two of the H iterations that the program prints in determining the R tQ network in Fig I. 
R=1.07kCl and RS=tUW kSl are the final mines. 



Circuit converts 7-segment code to decimal 



P Mazumder 

University of Illinois, Urbana, IL 

The 3-IC circuit in Fig 1 accepts 7-segment input code 
representing either of two commonly used numeric- 
character sets, converts it to an intermediate BCD 
code, and converts that to i 



code, defined in Fig 2, can be active high or active low. 
In addition, for convenience sake you can substitute 
exclusive-NOR (XNOR) gates for the exclusive-OR 
(XOR) gates shown. 

Table 1 shows the active-high output (Pj) that corre- 
sponds to each 7-segment character in the two sets, 
along with the intermediate BCD code (ABCD) that 
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